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Abstract 



We investigate numerically the transverse and longitudinal correlation lengths of the 
three-dimensional 0(4) model as a function of the external field H . From our data 
we calculate the scaling function of the transverse correlation length, and that of the 
longitudinal correlation length for T > Tc. We show that the scaling functions do not 
only describe the critical behaviours of the correlation lengths but encompass as well 
the predicted Goldstone effects, in particular the if~^/^-dependence of the transverse 
correlation length for T < T^.. In addition, we determine the critical exponent 
5 = 4.824(9) and several critical amplitudes from which we derive the universal 
amplitude ratios = 1.084(18), Qc = 0.431(9), = 4.91(8), = 1.265(24) 
and [/| = 1.99(1). The last result supports a relation between the longitudinal and 
transverse correlation functions, which was conjectured to hold below but seems 
to be valid also at Tc. 
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1 Introduction 



In 0{N) spin models with > 1 two types of correlation lengths appear. They 
govern the exponential decay of the correlation functions of the transverse and longi- 
tudinal spin components, defined relative to the external field H. Like in the case of 
the magnetization and the susceptibilities, the behaviour of the correlation lengths 
in the critical region is described by asymptotic scaling functions, critical exponents 
and amplitudes, which characterise the underlying universality class. These quanti- 
ties are of general interest. In addition, there are predictions PIE] for the correlation 
lengths, which are related to the presence of massless Goldstone modes jHl E] and are 
still untested. The measurement of the correlation lengths as functions of the field 
H enables us to verify these predictions and to determine the critical parameters 
and scaling functions. We have chosen the three-dimensional 0(4) model to carry 
out this program for the following reasons. Firstly, this model is of importance for 
quantum chromodynamics (QCD) with two degenerate light-quark flavours at finite 
temperature, because it is believed jSj-jni to belong to the same universality class 
as QCD at its chiral transition in the continuum limit. Secondly, the if-dependent 
model has already been investigated in some detail by Monte Carlo methods, in 
particular the critical behaviours of the magnetization jTOl E] and susceptibilities. 
In addition the corresponding Goldstone-mode effects have been verified jll . 

The specific model which we study here is the standard 0(4)-invariant nonlinear 
cr-model, which is defined by 

(3n = -J J2 ^--^y - H^-^ (1) 

<x,y> X 

where x and y are nearest-neighbour sites on a three-dimensional hypercubic lattice, 
and (f)g is a four-component unit vector at site x. It is convenient to decompose the 
spin vector (p^ into longitudinal (parallel to the magnetic field H) and transverse 
components 

<fe = 4eW + , with CH = H/H . (2) 
The order parameter of the system, the magnetization M, is then the expectation 
value of the lattice average 0" of the longitudinal spin components 

^ = = (0")- (3) 

X 

Here V = and L is the number of lattice points per direction. There are two 
types of susceptibilities. The longitudinal susceptibility is the usual derivative of the 
magnetization, whereas the transverse susceptibility corresponds to the fluctuation 
per component of the lattice average 0^ of the transverse spin components 

= H = V({*'=>-Af=), (4) 
X. ^^(^->^f. (5) 
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The expectation value (</>"'") is of course zero. The connected two-point correlation 
functions of the longitudinal and transverse spins are defined by 

Gl{S) = {44) -M\ (6) 

Gt{x) = 1(4- 4)- (7) 

They are related to the susceptibilities by 

Xl,t = J^Gz^.tI^) . (8) 

For all temperatures T (the coupling J acts here as inverse temperature, i. e. J = 
1/T) and fields H, except on the coexistence line H = 0,T < Tc and at the critical 
point, the large distance behaviour of these correlation functions is determined by 
the respective exponential correlation lengths ^l^t 

Gl,t{x) r~^exp{-\x\/^L,T) . (9) 

On the coexistence line, where the correlation functions decay according to a power 
law, it is still possible to define a transverse correlation length |T2] from the so-called 
stiffness constant. We do however not consider this option here. 



1.1 Goldstone-Mode Effects 

The spontaneous breaking of the rotational symmetry for temperatures below the 
critical point gives rise to the so-called spin waves: slowly varying long-wavelength 
spin configurations, whose energies may be arbitrarily close to the ground-state en- 
ergy - the massless Goldstone modes ^HI- For H ^ the magnetization M below 
Tc attains a finite value, the spontaneous magnetization M{T, 0) (here and in the 
following we assume always that H > 0, so that M{T, 0) > 0). The transverse sus- 
ceptibility xt = M/ H, which is directly related to the fluctuation of the Goldstone 
modes, diverges therefore as when /J ^ for all T < Tc. It is non-trivial that 
also the longitudinal susceptibility xl is diverging on the coexistence curve. The 
predicted divergence in three dimensions is [H] 

XL{T<Tc,H)r^H-'/^ . (10) 

From a phenomenological spin wave analysis and the behaviour of the transverse 
susceptibility xt Fisher and Privman pP arrived at the conclusion that the bulk cor- 
relation length, which they identify with the transverse correlation length, diverges 
also when if as 

^t{T<Tc,H)^H-^'^ . (11) 

In fact, this connection between the behaviours of the correlation length and the 
respective susceptibility, namely 

el ~ XT, (12) 



2 



is what one usually expects ^]. For instance, when masses are defined from the 
corresponding inverse susceptibilities, 



m' = , (13) 

one would anticipate that m ~ C,^^. In QCD, the transverse mass corresponds [7j 
to the pion mass, rriT = m.„, and the longitudinal one to the sigma mass, = m„. 
According to Ref. P, the relation equivalent to Eq. (fT^ . C,1 ^ Xl, does not hold 
for the longitudinal fiuctuations at long wavelengths, because they are driven by 
the transverse fiuctuations and the bulk correlation length also sets the scale of the 
decay of Gl{x). This statement is substantiated by the well-known relation jU I2I 
between the longitudinal and transverse correlation functions at zero field, T < 
and large distances |x| 



Gt{x) 
M 



1 2 



(14) 



where in our case = 4. One expects that the relation is still valid for small 
non-zero fields H near the phase boundary in the region of exponential decay. That 
implies a factor of 2 between the correlation lengths. 

The rest of the paper is organized as follows. First we discuss the critical be- 
haviours of the observables and the universal scaling functions, which we want to 
calculate. In Section 3 we describe some details of our simulations and the way we 
determine the correlation lengths. Section 4 serves to find the critical amplitudes, 
which are needed for the normalizations and the universal ratios. In the following 
Section 5 we discuss the scaling functions which we obtain from our data. The 
Goldstone effect on the transverse correlation length is demonstrated in Section 6. 
Subsequently we investigate the if-dependence of the correlation lengths in the high 
temperature phase. We close with a summary and the conclusions. 



2 Critical Behaviour and Scaling Functions 

In the thermodynamic limit {V oo) the observables show power law behaviour 
close to Tc. It is described by critical amplitudes and exponents of the reduced 
temperature t = {T — Tc)/Tc. The scaling laws at if = are for 

the magnetization 

M = B{-tf for t < , (15) 
the longitudinal susceptibility 

XL = C+r^ for t > , (16) 
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and since for H = 0,t > the correlation lengths coincide = = ^ {like the 
susceptibilities) 

^ = ^+1-" for t > . (17) 
On the critical line T = Tc or t = we have for H > the scaling laws 

M = B^H^'^ or H = D^M^ , (18) 

XL = C^H^^^'^ with = By 5 , (19) 
and for the correlation lengths 

a,r = Cl,tH'''^ , = vlfi^ . (20) 

We assume the following hyperscaling relations among the critical exponents to be 
valid 

7 = /3(5-l), rfz/ = /3(l + 5), 2-^ = ^|v. (21) 

As a consequence only two critical exponents are independent. Because of the hyper- 
scaling relations and the already implicitly assumed equality of the critical exponents 
above and below one can construct a multitude of universal amplitude ratios [12] 
(see also the discussion in Ref. JE])- The following list of ratios contains those which 
we will determine here 

= C+D,B'-' , Ul = Ct/Cl , (22) 

= B\er/c+ , Q2'^ = {CL,T/er^''c^/c'' . (23) 



The critical behaviour of the magnetization in the vicinity of is more generally 
described by the magnetic equation of state. In its Widom-Griffiths form it is given 

by mm 

y = fix) , (24) 

where 

y = h/M^ , X = i/M'^"^ . (25) 

The variables i and h are the normalized reduced temperature i = ITc/Tq and mag- 
netic field h = H/Hq, which are chosen such as to fulfill the standard normalization 
conditions 

/(O) = 1 , /(-I) = , (26) 

which imply 

M{t = 0) = /i^/^ and Ho = D,, (27) 

M{h = 0) = (-t )^ and To = B-^/^T, . (28) 

Possible dependencies on irrelevant scaling fields and exponents are however not 
taken into account in Eq. ()24|) . the function f{x) is universal. Another way to 
express the dependence of the magnetization on i and h is 

M = h}'^fG{tlh}/^^) , (29) 
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where fc is again a universal scaling function. The two forms (j^^ and (j^^ are of 
course equivalent. The function fciz) and its argument z are related to x and y by 

fG = y-'/', z^i/h'l^' = xy-''^' . (30) 

Correspondingly the normalization conditions translate into 

/g(0) = 1 , and fa{z)^-^J-zf . (31) 

Since the susceptibility xl is the derivative of M with respect to H we obtain from 
Eq. 

with 

/x(^) = \ [fciz) - ^/a^)) • (33) 

For if — >■ at fixed t > 0, that is for 2; — * oo, the leading asymptotic term of is 
determined by Eq. (fTBj) 



For z ^ 00 the leading terms of fa and are identical, because for T > Tc and 
small magnetic field M is proportional to if. 

Like for M and xl the dependence of the correlation lengths on i and h in the 
critical region and the thermodynamic limit is given in terms of scaling functions 

gl''^{z) by 

= h~''^gf{z) . (35) 

These functions are universal except for a normalization factor. On the critical line 
t = or 2; = we find from (j^Ull 

gfi^) = Cl,tD:''^ = ilABy'^ , (36) 

and from ()17p the asymptotic behaviour at z ^ 00 (which is the same for both 
correlation lengths), 

^?''^(^) .^00 tB^'^z-^ . (37) 
Indeed, the ratios of the amplitude for z ^ 00 in and the g^'^ (0) are universal 



(3J 



whereas the 5f|''^(0) itself are not. In analogy to the Ising case discussed in Ref. [T 
we therefore define two universal scaling functions by 



gf{z) = gf{z)/gf{Q) . (39) 
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3 Numerical Details 



All our simulations were done on three-dimensional lattices with periodic boundary 
conditions and linear extensions L = 48, 72, 96 and 120. As in Ref. [H^ we have 
used the Wolff single cluster algorithm jTH!, which was modified to include a non- 
zero magnetic field [20]. In order to reduce the integrated autocorrelation time Tint 
we performed 100-3000 cluster updates between two measurements for H > 0, such 
that we obtained Tint ^ 3 for L = 48 and < 5 for L > 72. For H = 0,T > Tc 
and large lattices the number of cluster updates was increased up to 10000. In 
general we made 20000 measurements for each fixed H and J. We use the value for 
Jc = = 0.93590 from Ref. ^T]. The coupling constant region which we have 
explored was 0.9 < J < 1.2, the magnetic field was varied from H = to H = 0.01. 



3.1 Measurement of the Correlation Lengths 

Instead of using correlation functions of the individual spins it is more favourable 
to consider spin averages over planes and their respective correlation functions. For 
example, the spin average over the {x, ?/)-plane at position z is defined by 

^^ = J2T.^^^- (40) 

x,y 

The spin averages have again longitudinal and transverse components 

= Slen + S^ , (41) 

The expectation values of the components are independent of z and equal to those 
of the respective lattice averages 

(^^) = = M , and {S^) = {$^) = 0. (42) 

Correspondingly, we define now plane-correlation functions Gl^t{z) by 

Gl{z) = (^{SlsD - M') (43) 

Gt{z) = ^{S^S^). (44) 

Here, z is the distance between the two planes. Instead of choosing the 2;-direction 
as normal to the plane one can as well take the x- or ?/-directions. Accordingly, we 
enhance the accuracy of the correlation function data by averaging over all three 
directions and all possible translations. The correlators are symmetric and periodic 
functions of the distance r between the planes, the factor on the right-hand sides 
of (PHj) and ensures the relation 

L-l 

XL,T = J2^LAr) ■ (45) 

r=0 
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Figure 1: The effective transverse correlation length ^^■^■^ at J = 0.95 and H = 0.001 
for L = 48 (crosses), 72 (diamonds), 96 (squares) and 120 (stars) as a function of r. 
The lines show the corresponding r-ranges where fits to Eq. ()4(i|l were performed, 
their heights the resulting fit values. The data for L = 72, 96 and 120 have been 
shifted upwards by 1 , 2 and 3 for better visibility. 



Like the point-correlation functions in Eq. Q the plane correlators Gl^t{t) decay 
exponentially. In order to obtain the corresponding exponential correlation length 
^ from G{t) we therefore make the ansatz 



G{r) = A [ exp(-r/0 + exp(-(L - r)/0 



(46) 



and then try to fit the data for the correlation functions in an appropriate r-range. 
The ansatz ()46|1 implies of course, that there are no additional excitations con- 
tributing to G{t). Inspired by the experiences reported in the Ising case [221 
we proceed in the following way. First we calculate an effective correlation length 
^^^■^{t) from ()46|) using only the correlators at r and t+1. For t <^ L this correlation 
length is approximately given by 



-1 



ln(G'(r+ 1)/G(r)) 



(47) 



As an example, we show in Fig.[T]a typical result for ^^^^ {t) at J = 0.95, H = 0.001, 
that is in the low temperature region, from the transverse correlation functions on 
lattices with L = 48, 72, 96 and 120. With increasing r also ^^-^-^ increases and even- 
tually reaches a plateau inside its error bars. Most likely, the lower values of ^'^■^■^ at 
small r are due to higher excitations. At large distances r the resulting ^'^■^•^ start 
to fluctuate when the relative errors of the data become too large. The two limits 
define an intermediate r-range where a global fit with the ansatz can be used 
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0.93 


11.789(36) 


11.813(32) 


6-26 


7.594(57) 


7.647(60) 


8-19 


0.92 


8.563(32) 


8.592(26) 


3-17 


7.705(54) 


7.720(60) 


4-15 



Table 1: Correlation length estimates for different J and H = 0.001: ^'^^^{to) is 
the value at the lower end of the plateau region [r^, rf,], ^ the value from fits to the 
ansatz in the whole plateau region. 



to estimate the exponential correlation length ^. It is clear, that due to the strong 
correlations between the values G{t) for different r, the result from a fit in the whole 
range of the plateau must be essentially the same as that of one (low) r-value in the 
plateau. That is indeed the case: the two estimates for ^ agree within the error bars 
and the errors themselves are of about the same size, or somewhat larger for the 
plateau fits. Using the latter makes the result less dependent on the special choice of 
r, where ^^^^ is evaluated, in case the ^'^■^■^-values are slightly fluctuating. In order 
to give an impression of the different ^-estimates, their errors and the respective 
r-ranges we show in Table Q some representative results for T below, at and above 
Tc and H = 0.001 on a lattice with L = 72 (at J = 0.92 with L = 48). 

With our method we found reliable values for the exponential transverse cor- 
relation length C,T at all T. In the case of the longitudinal correlation length the 
same is true in the high temperature phase T > Tc. However, at Tc and slightly 
below the critical point at if > there are already indications for the presence 
of higher states. Deeper in the low temperature phase it becomes impossible to 
estimate the longitudinal correlation length: no plateau appears, there is a dense 
spectrum of states contributing to the longitudinal correlator. 



4 Critical Amplitudes and Universal Ratios 

In Ref. the spontaneous magnetization M(T, 0) was calculated at several tem- 
peratures on the coexistence line, taking advantage of the Goldstone effect on M. 
From these results the critical amplitude B of Eq. (|T5|) and the normalization Tq 
could be determined to B = 0.9916(5) and Tq = 1.093(2) using f3 = 0.380. We use 
these values in the following. In the same paper the amplitude B'^ and the critical 
exponent 6 were determined from fits to the critical behaviour of M on the critical 
line. Since we have done new simulations on the critical line with higher statis- 
tics and at more if-values on the larger lattices than in ^l] we could improve the 
analysis with the new magnetization data. 
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Figure 2: The magnetization M at as a function of the magnetic field H for 
L = 48, 72 and 96. The solid line shows the fit (jl8|l . 

4.1 Results from the Critical Line 

On the critical isotherm, that is at Tc, we have measured the magnetization for 
H > on lattices with L = 48, 72 and 96. The results are shown in Fig. |21 We 
observe essentially only for L = 48 a noticeable finite size dependence close to H = 0. 
Like in [llj we have made various fits to the data from the largest lattices in the 
if-interval [0.00075, 0.00425] with the simple ansatz 

M = B^H^/^ . (48) 

The resulting values for the amplitude and the exponent are (with /Nf 1.7) 

B" = 0.721(2) and S = 4.824(9) , (49) 

and lead to Dc = Hq = 4.845(66). The solid curve in Fig. |21 shows this fit. It also 
represents the data at higher magnetic fields very well. From these additional data 
a marginal negative correction-to-scaling term may be inferred, which is however 
irrelevant at low H. Compared to the corresponding fit in jTT], where S = 4.86(1), 
we note that our new result for S is closer to the value 4.789(6) from Ref. |23j. The 
latter value was however deduced from calculations at zero magnetic field. All the 
other critical exponents which were needed in our calculations have been derived 
from the hyperscaling relations (PT|) and (3 = 0.380, 6 = 4.824 . We use therefore 

u = 0.7377 , 7 = 1.4531 , i^c = 0.4024 . (50) 

The data for the transverse and longitudinal correlation lengths on the critical line 
are shown in Fig. 01 (a) and (b). For the transverse correlation length C,t we find only 
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Figure 3: The correlation lengths (a) and C,l (b) at for L = 48, 72 and 96 as a 
function of H. The lines are the fits ()5H) with C,t l froni (|52|l and ct,l = 0. 

finite size effects very close to H = 0. The longitudinal correlation length C,l exhibits 
fluctuations and a systematic deviation to higher ^^.-values, when the magnetic field 
is decreasing. The smaller the lattice, the earlier this behaviour sets in. In order 
to determine the amplitudes we have fitted our results to the following form with 
u = 0.8 

^T,L = eT,LH-''^a + CT,LH-''^) . (51) 
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In both cases we have taken various subsets of the data for the fits in the iY-interval 
[0.0005,0.005]. It turned out that the correction terms are zero within their error 
bars and that fits with ct.l = work just as well. The X^/^f varies between 0.5 
and 0.8 . We find for the amplitudes 

Ct = 0.838(1) , a = 0.421(2) . (52) 

As can be seen from Fig. El (a) and (b), the corresponding fits describe the data 
also in the higher if-range up to 0.01 . The ratio of the two correlation lengths is 
remarkably close to 2 

Ul = CtICl = 1-99(1) , (53) 

that is, the expectation for the ratio from Eq. ()14|) close to the coexistence line is 
fulfilled also at T^. A similar observation has been made for the three-dimensional 
0(2) model |2lj. 



4.2 Results for T > and H = 

Our next aim is the determination of the critical amplitudes and To this 
end we have evaluated the data of the susceptibility and of the correlation length 
above the critical temperature for H = 0. The data points for x = Xl = Xt are 
plotted in Fig. |3] as a function of t. We have made various fits with the ansatz 

X = C+t-^[l + C+n , (54) 
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Figure 5: The correlation length ^ above T^. at H = for L = 48, 72, 96 and 120 
versus t. The solid line shows the fit from Eq. ()57|1 . 



in the J-interval [0.91,0.932], which corresponds to 0.00415 < t < 0.0285. All fits, 
including a correction term ~ or not, are compatible with the critical amplitude 
result 

C+ = 0.231(2) . (55) 

For the correction amplitude we find = 0.14(7), and x^/^f is in the range 
[1.3, 1.8] . As can be seen in Fig. IH where only the leading term is plotted, the 
correction contribution is marginal. 

The correlation length data for = above the critical temperature are shown 
in Fig. El Like for the susceptibility we have made an ansatz including a correction- 
to-scaling term 

^ = ^+t-^[l+^+r^] . (56) 

Different fits in the same J-interval [0.91, 0.932], which was already used in the case 
of the susceptibility, lead to the result 

e = 0.466(2)t-'^[l - 0.12(4)t'^^] , (57) 

with a /Nf between 0.8 and 1.2 . The critical amplitude of the correlation length 
is therefore 

^+ = 0.466(2) . (58) 
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4.3 Universal Amplitude Ratios 

The critical amplitudes B'^, ^£5., C'^, which we have just determined, and B 
from pTj enable us to fix the necessary normalizations of the scaling functions and 
in addition we can calculate the universal amplitude ratios defined in Eqs. (12 2j) and 
f)2Hj] ■ The ratio f/| has already been discussed. We obtain for the other ratios 

= 1.084(18) , Qc = 0.431(9) , = 4.91(8) , = 1.265(24) . (59) 

For the ratio Ry. several alternative values exist: from the parametrization of the 
equation of state in ^l] one gets 1.126(9), from another parametrization in [HI 
1.12(11), from the 1/iV-expansion of Oku and Okabe 1.098, and from the e- 
expansion 1.239 ^2]- Our finding for Qc is well in accord with the result 0.44(2) of 
j9j. We could not find competing values for the ratios Q2' t^e literature. 



5 The Scaling Functions 

In Ref. PI] the equation of state in the Widom- Griffiths form was parametrized by 
a combination of a small-x (low temperature) form Xs{y), which was inspired by the 
approximation of Wallace and Zia Jl] close to the coexistence line (x = —1; y = 0) 

xs{y) + 1 = {ci + d^)y + Z2y^'^ + d^y^'^ , (60) 

and a large-x (high temperature) form xi{y) derived from Griffiths's analyticity 
condition [TH] 

xi{y) = ay'/'' + hy^^'^^^'^ . (61) 
The two parts are interpolated smoothly by the ansatz 

xiy) = Xs{y) „ + xi{y) , (62) 

from which the total scaling function is obtained. Since we have changed the val- 
ues of Hq and of 5 in comparison to Ref. [Hj, we have redone the corresponding 
parametrization fits. With the new data we find 

ci + = 0.359(10) , C2 = 0.666(6) , (63) 

and = 1 — (ci + 6/3 + C2) is fixed by the normalization condition ?/(0) = 1. The 
new values for a and h are 

a = 1.071(4) , b = -0.866(38) . (64) 

The interpolation parameters yo = 10.0 and p = 3 have been retained unchanged. 
The modified values of a and 7 result in a new estimate of R^ 

R^ = a? = 1.105(4) , (65) 
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hne shows our parametrization. The numbers refer to the different J = 1/T- values 
of the data, the dotted hues cross at the normahzation point /g(0) = 1. 
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Figure 7: The scahng function fx{z) = xlHqK^'^^^ versus z = The sohd 

hne is from our parametrization of the equation of state. The numbers refer to the 
J = 1/T-values of the data. The normahzation is /;^(0) = 1/6. 
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in nice agreement with the value we obtained in directly from the amplitudes. 

In Fig. iniwe present the results for fdz) obtained from the magnetization data 
in the coupling range 0.90 < J < 1.0. At first sight these data are scaling well. Since 
the data from the coupling J = 1.2 showed already visible deviations from scaling 
behaviour we have discarded them. Similarly, limitations were found for the scaling 
if-regions. The data for fciz) were scaling in the low temperature region only up 
to H = 0.005 and in the high temperature region for J = 0.93 and J = 0.92 up to 
= 0.05, for J = 0.91 and J = 0.90 the scaling if-range extended only to if = 0.01. 
Only data within these parameter ranges are shown in Fig. IHl. In the same figure we 
show the function fciz), which one obtains with Eq. ()30|1 from our parametrization 
for x{y). It obviously describes the data quite well. The scaling function f^iz) 
of the susceptibility is connected via Eq. to the scaling function fciz) of the 
magnetization. It can therefore be calculated from the parametrization as well. On 
the other hand we have direct data for xl from our simulations which allows for 
another check of the scaling hypothesis and a comparison to the parametrization. 
In Fig. [3 we show the respective data for the same J- values as in Fig. IHl The data 
are not as accurate as those for the magnetization. We observe explicit scaling in 
the high temperature phase {z > 0), but with decreasing temperature, or larger —z, 
the data are spreading and no longer scaling. In particular the data at the coupling 
J = 1.0 seem to be already outside the critical region. For J = 0.975 and J = 0.95 
it is possible that at very small H (large —z) finite size effects are responsible for 
that behaviour. The other data are coinciding in essence with our parametrization 
for z > —3. 



5.1 The Scaling Functions of the Correlation Lengths 

In Fig. IHlwe show the data for the normalized scaling function gj{z) = ^rh'"" / gj [0) . 
The determination of the transverse correlation length with the method described 
in Section 13.11 was always possible. In each case a plateau in ^^■^■^ could be found. 
The data shown in Fig. |H1 correspond to the ones presented in Fig. [7| for the scaling 
function of the susceptibility, apart from the points for J = 1.0 , which were clearly 
outside the scaling region. The normalization was calculated from Eq. ()36p and is 
gj{0) = 0.444(3). As in the case of the magnetization and the susceptibility we 
observe scaling of the correlation length data in a large z-range. For 2; > we see an 
early approach to the asymptotic form calculated from Eq. (jTTj) . The shape of gj{z) 
is very similar to that of faiz)- It is not by accident that we find this behaviour 
but it is due to the relation (fT^ between the transverse correlation length and its 
susceptibility, which we discussed in Section IlD . We consider therefore the ratio of 
and xt and find 

XT IG 

the same relation as in the Ising case except that here is replaced by fc, 
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Figure 8: The scaling function of the transverse correlation length gj{z) = 
irh"" / gj {Q) . The solid line is the asymptotic form calculated from Eq. pTjl . The 
numbers refer to the different J = 1/T- values of the data. 



which plays the corresponding part for the transverse correlation length. In the 
Ising case we found the asymptotic behaviour 

^ ~ i±z)-'"' . (67) 

JX 

It is easy to show that {gJY/ fa for z +oo is proportional to z~'^'^. If for 
the transverse correlation length behaves as ~ H~^/'^ - as expected because of 
the Goldstone effect - then for z —>■ — oo the ratio [gj^/fc must be proportional 
to {—z)~'^'^, as in the Ising case. The condition for the asymptotic behaviour of the 
ratio in the broken phase and the if~^/^-behaviour of are thus equivalent. 

Assuming the corresponding asymptotic behaviour for {g^Y / fx would imply that 
~ H^^/'^ for T < Tc, in contradiction to the expectation of a factor of 2 between 
the two correlation lengths. The situation in the low temperature region is indeed 
difficult and as it seems we cannot clarify numerically the status of there. In Fig.lHl 
we show our results for the normalized scaling function g^{z) = C.ih"" / g^{0), where 
g^{0) = 0.223(2) was determined from Eq. (jHUj) . In the high temperature region 
the longitudinal correlation length could be calculated with sufficient accuracy, the 
finite size effects were negligible. Approaching the critical temperature {z = 0) from 
above, the data become more and more noisy and below Tc we are essentially unable 
to determine and g^{z): if one can find a plateau in ^^■^■1' at all, then the results 
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Figure 9: The scaling function of the longitudinal correlation length g^{z) = 
^ih^" /g^{0)- The solid line is obtained from the asymptotic form ()H7|1 . The numbers 
refer to the different J = 1/T- values. The dashed line indicates the peak position 
of xl, the symbol notation is as in Fig. |H1 



for different lattice sizes are often very different. In addition scaling is completely 
lost already for small negative values of z. We have therefore only included the data 
for z > —2 in Fig. El to show this behaviour. As compared to the transverse scaling 
function in the high temperature region we have a different functional form, which 
is similar to that of the scaling function of the longitudinal susceptiblity in Fig. [7| 
Like in the Ising case, both functions have a peak at about the same z-value, which 
is at Zp ~ 1.335 in f^, in accord with the value 1.33(5) found in Ref. [Hj. As in Fig. 
IHlwe find an early approach to the asymptotic form given by Eq. P7jl . 



6 The //-Dependence of at Fixed T < 

Goldstone-mode effects are expected for all temperatures below Tc, not only in the 
critical region, but also there. We have discussed already in the last section, that 
the shape of the scaling function for the transverse correlation length is explicitly 
determined in its asymptotic part for z — oo by the predicted dependency of 
on for if — >■ 0. We have tested this prediction by calculating at three 

fixed temperatures below Tc. The data for are plotted versus 
for J = 0.95,0.975 and J = 1.0 . Actually, the same data points from the larger 
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Figure 10: The transverse correlation length as a function of H^^^^ for J = 
0.95, 0.975 and 1.0 . The lowest values are those of J = 0.95, the highest belong 
to J = 1.0. For better visibility we have lifted the results for J = 0.975(1.0) by 
5.0 (10.0). The dashed lines represent the straight hne fits to the data. 



lattices, apart from those for J = 1.0, had already been included in our scaling plot. 
We see from Fig. El that the data always follow a straight line 

= xo{J)H-'/^ + xi{J) , (68) 

where the slope is slightly increasing with increasing J: xq = 0.481(2), 0.577(1) and 
0.637(2) for the three J-values. The temperature dependence of Xq can be derived 
from the scaling function by comparing its leading asymptotic behaviour to the 
if-dependence of ^t- H 

gji^) = 9<ii-zy + . . . , (69) 

with gg a constant, then the power q must be g = P6/2 — u to lead to the required 
if-dependence for ^t- That implies 

xo = 9,{~i)m',/' , (70) 

and since q = 0.1789 is positive, we have an increase with increasing J. We can 
check the J-dependence of the slopes xq by taking ratios 

^o(-^i) f ^^^^ 



Xq{J2) \-t 

For xo(0.975)/xo(0.95) = 1.200(5) our formula predicts 1.195, that is rather close to 
the fit value. For xo(1.0)/a;o(0.975) = 1.104(4) we get 1.088, a slightly lower result. 
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On the other hand, the data for J = 1.0 were just outside the scahng region and 
somewhat higher than the scahng function. The small deviation to a higher ratio is 
therefore what one expects. 



7 The //-Dependence of ^t,l at Fixed T > Tc 

At H = the correlation lengths in the symmetric phase (T > Tc) are equal. We 
have shown the t-dependence in Fig. El With increasing H, however, they differ 
and is always bigger than ^l- In Fig- ^2 we show the two correlation lengths as 
a function of H for three different fixed couplings or temperatures: J = 0.92,0.91, 
and 0.90. The lattice sizes were L = 48 and 72. Essentially no finite size effect was 
found. As expected, the correlation lengths decrease with increasing field from their 
value at H = 0. The reduction itself is diminishing with increasing temperatures 
(decreasing couplings J = 1/T), the curves become flatter and the difference between 
the two correlation lengths disappears. This behaviour can be understood from the 
asymptotic expansions of the scaling functions of the correlation lengths. In the 
symmetric phase the correlation lengths are even functions of h, so that C,{h) = ^{—h) 
(we omit the indices T and L for this discussion). The scaling function g^{z) must 
then have an asymptotic expansion for z ^ oo {h ^ 0,t fixed) of the form 

oo 

g^{z) = z-^Y.9nZ-'^\ (72) 

ra=0 

where the gn are constants, go = ^^B'^^^, and A = pS. This translates into the 
following /i-expansion for the correlation length at fixed t and small h 

oo 

^ = {t)-''J29n-{i)-'''''h'- . (73) 

n=0 

For very small magnetic fields h we have 

e = goii)-' + g,{iy'-'''h' + ... , (74) 

and with increasing i the /i-dependence will disappear with a factor ~ it)~^^ com- 
pared to the value at h = 0. 

In the neighbourhood of the critical point, where \z\ is small (but not h), one 
can expand g(^{z) in powers of z and obtains 

e = h-^^g^{z) = h-^^ [gl + glz + g^z^ + . . .) , (75) 

with g^ = ^^[B'^y/^ . That leads to the following /i-dependence of ^ 

i = h-^^ {g^, + glih~'/^ + gm'h-'/^ + ...), (76) 



19 



10 



5 



1 

- 1 






1 ' 

^ T 7, 
± f 1-1 




' 1 

- 


' ▼ — ^^-▼=^ 




— ▼ 






- 


:| 






T 

V 


▼ 








■ 


■ 


— _ 

■ 


' S7 


i 1 ... 1 



0.001 0.005 0.010 

Figure 11: The correlation lengths C,t (filled symbols) and (empty symbols) as 
a function of H. The couplings are J = 0.92 (circles), J = 0.91 (triangles) and 
J = 0.90 (squares). In order to disentangle the different curves we have lifted the 
results for J = 0.92 (0.91) by 5.0 (2.5). The data are connected by straight lines to 
guide the eye. 

when the temperature is near to the critical point. In Fig.^Jwe show as an example 
the correlation lengths for J = 0.93, close to Jc = 0.93590 . We have made fits in 
the if-interval [0.0006,0.005], using only the first two terms in Eq. (|76|l . Obviously 
these simple fits describe the data quite well in a large range. 



8 Summary and Conclusions 

In this paper we have investigated the correlation lengths of the three-dimensional 
0(4) model with Monte Carlo simulations. One of the main objectives was to find 
the functional dependence of the correlation lengths on the external field at fixed 
temperature, a second to understand the interplay of the critical behaviour and the 
effects induced by the massless Goldstone-modes. A main result of our work is the 
calculation of the scaling function of the transverse correlation length and that of the 
longitudinal correlation length for T > Tc. In the low temperature phase we were. 
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Figure 12: The correlation lengths (upper set of data) and (lower set) as a 
function of H for the coupling J = 0.93. The lines are fits to Eq. (|76|l in the interval 
[0.0006,0.005]. 

however, unable to reliably estimate ^l- The reason for that is the rich spectrum of 
states which contribute to the longitudinal correlators below the critical point. In 
the simpler Ising model where only one correlation length exists, the influence 
of higher states [SE] complicates the determination of the correlation length below 
Tc, but in contrast to the 0(4) model it is still possible there. 

As we could show, the scaling functions do not only describe the critical be- 
haviours of the correlation lengths, but encompass as well the predicted Goldstone 
effects. Indeed, in consequence of the relation ~ xt and the equation xt = M/ H 
for the fluctuation of the transverse spin components, the prediction ~ H^^^^ 
emerges below Tc and leads to a similar functional form of the scaling functions of 

and that of the magnetization. In the high temperature phase a similar corre- 
spondence exists between the scaling functions of and that of the longitudinal 
susceptibility. There both functions have a peak at the same position, a behaviour 
which was found also in the Ising model [T3] . 

In preparation of the scaling functions we have determined several critical ampli- 
tudes of the correlation lengths: ^£ and and in addition also B'^ and C~^. As 
a byproduct we have found the critical exponent 6 = 4.824(9). The critical ampli- 
tudes allowed us to directly determine the universal amplitude ratios R^, Qc, Q2' 
and f/| = 1.99(1) . The latter result is most remarkable, because it confirms a 
relation between the longitudinal and transverse correlation functions, which was 
conjectured to hold below Tc, but seems to be valid also at Tc. 
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